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We design an adiabatic quantum algorithm for the counting problem, i.e., approximating the 
proportion, a, of the marked items in a given database. As the quantum system undergoes a 
designed cyclic adiabatic evolution, it acquires a Berry phase 2-7ra. By estimating the Berry phase, 
, we can approximate a, and solve the problem. For an error bound e, the algorithm can solve the 

problem with cost of order (-|) 3 ^ 2 , which is not as good as the optimal algorithm in the quantum 
circuit model, but better than the classical random algorithm. Moreover, since the Berry phase is 
£Nj , a purely geometric feature, the result may be robust to decoherence and resilient to certain noise. 
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I. INTRODUCTION 



Quantum algorithms can solve certain problems significantly faster than the known classical algorithms. However, 
one main obstacle to realize the faster quantum algorithms is the decoherence induced by the coupling environment. 
To overcome the obstacle, a novel quantum computation model, quantum adiabatic computation, is a promising 
candidate [J [H, H, 0, H, H, 0, § • The model is believed to enjoy inherent robustness against the impact of decoherence 
in [1, H, @|. Even though not all researchers share this view @, adiabatic quantum computation still attracts 
considerable attention. It has been proved to be polynomially equivalent to the quantum circuit model [3, 0|. For 
instance, in the adiabatic model searching an unordered database requires time of the same order of magnitude 
, as Grover's algorithm 0, Q, but few adiabatic algorithms are known that have performance similar to that of the 
■ corresponding algorithm in the quantum circuit model. In this paper, we show an adiabatic algorithm for the counting 
t^- \ problem. The task of the counting problem is to approximate the proportion of marked items in an TV-item database, 
which is denoted by a. In classical computation, the counting problem needs O(N) evaluations in the worst case 
setting, and 0{{\) 2 ) in the randomized setting, for error e. There exists a quantum algorithm solving the problem 
in 0{\) evaluations [1(3]. Moreover, Nayak and Wu [ll| showed that the quantum algorithm given in [1(3] is optimal 
in the quantum circuit model. The countin g pr oblem is also central for many continuous problems, such as high 
dimensional integration, path integration 1^, IH [H| and eigenvalue approximation pj| . 

Our algorithm for the counting problem is based on the Berry phase acquired in the adiabatic evolution. When a 
quantum system undergoes a cyclic adiabatic evolution, it acquires a geometric phase, which is known as the Berry 
phase [16]. In the algorithm, we design an adiabatic evolution such that the resulting Berry phases encode the solution 
of the problem. Then we can solve the problem by estimating the Berry phase after the adiabatic evolution. Since 
Berry phases are global phases and cannot be measured directly, we let two parts of a superposition undergo the 
same cyclic adiabatic evolution in different directions. After the adiabatic evolution, the dynamic phases of the two 
parts cancel out, and the Berry phases are ±27ra. Then, by estimating the relative phase between the two parts of 
the superposition, we can estimate a. In the static model, the algorithm has a runtime of order (^) 3 ^ 2 , which beats 
the optimal classical algorithm in the randomized setting. Usually, in adiabatic algorithms it is the final state that 
encodes the solution of the problem, while in our algorithm it is the Berry phase. Since the Berry phase is a purely 
geometric feature, i.e., it only depends on the path of evolution, and is independent of details of how the evolution is 
executed, the result is resilient to certain small errors. 

We begin with some preliminaries about adiabatic algorithms and the Berry phase, which are helpful in presenting 
our main result. Consider a quantum system in a state \i/j(t)}, (0 < t < T), which evolves according to the Schrodinger 



'Electronic address: cz2165@columbia.edu 
tElectronic address: weich03@mails.tsinghua.edu.cn 
tElcctronic address: ap@cs.columbia.edu 



2 



equation 

l f t \m) = H(t)\m), (i) 

where H (t) is the Hamiltonian of the system at time t. If the system is initially in its ground state, and the Hamiltonian 
varies slowly enough, it will remain close to the ground state of H(t), at time t. Let \E (t)} be the ground state of 
the Hamiltonian, and Eo(t) be the corresponding eigenvalue. If H(T) = H(0), then \ip(T)) is close to \ip(0)}, with 
the exception of a global phase. The phase can be divided into two parts: the dynamic phase 

6 = - [ E {t)dt, (2) 
Jo 

and the geometric Berry phase 

1 = 1 1 ( E °\j t E ° )dt - (3) 

The Berry phase depends only on the path taken, not on how fast the path is traversed. Hence, if we design a cyclic 
path of Hamiltonians, the Berry phase is totally determined. 

The remainder of the present paper is organized as follows. In section UH we provide the basic adiabatic evolution 
used in the algorithm, which can encode the solution to a Berry phase of a quantum system. Then we give the 
adiabatic algorithm for the counting problem. In section IIII1 we will show the relationship between the accuracy of 
the algorithm and the time it used. In the Appendix, we show the detailed derivation of the difference of real relative 
phase and the Berry phase. 



II. THE ADIABATIC ALGORITHM FOR THE COUNTING PROBLEM 



In this section, we will show how to encode the solution to the Berry phases. We use a function / : {0, • • • , N— 1} — > 
{0, 1} to denote whether an item is marked, i.e., f(s) — 1, if the s-th item is marked; f(s) — 0, otherwise. For an 
error bound e, the algorithm has m = log(i) adiabatic evolutions. In each evolution, we use n + 1 qubits, where 
n = log TV. The quantum state can be divided into two systems: the control system, which has 1 qubit, and the 
computing system, which has n qubits. The computing system is in an iV-dimensional Hilbert space, whose basis 
states are denoted as \k), where k = 0, • • • , N — 1. We use an equal superposition of all basis states as the initial state 
in the computing system 

1 

^ k=0 

and use ^(|0) + |1)) as the initial state in the control system. Hence, the initial state of the whole system is 

|^0)) = -±=(|0) + |l))®hAo). (5) 
In each evolution, the initial Hamiltonian of the computing system is 

H = I-\tpo){H- (6) 
Typicallyin quantum algorithms, the fundamental oracle used is 

\x)\y)^\x)\y®f(x)). (7) 

In our algorithm, we modify the oracle to 

\x)\y)\z)->\x)\y(Bf(x))\z®(vf(x))). (8) 

With a 2-qubit auxiliary state -^(|0) ® ^ (|0) - |1)), using the oracle on the initial state, and then discarding 

the auxiliary qubits, the operation of the oracle can be written as 

1 N-l N-l 

7^ E l fc ) - 7^ E exp(-4/(fc))|fc). (9) 

v fc=0 v fe=0 
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Let j3 = 1 - a, M = aN, and 



|6)= ; J" \S), 

V s:/(«)=0 

ii> = -^ y |«>, 

s:/(s)=l 



(10) 



the initial state can be rewritten as 



l^o) = v^|6) + V5|i>, (ii) 

and the states after repeatedly using the oracle are 

m = v^|6) + H)\/S|i>, (12) 

for k = 1,2,3. 

In the adiabatic algorithm, we define four Hamiltonian oracles as 

# fe = /-hWfc|, (is) 

for k — 0, 1, 2, 3. Then consider a linear interpolation between the four oracles, that is 



H(t) = J2s k (t)H k , (14) 



fc=0 



where ^ Sk(t) = 1, for any < t < T. 
In the j-th evolution, we choose 



s = ^(1 + cosfy), si = ^ sin^j, 



-(1 - cos^),s 3 = --; 



(15) 



where is a function from [0, T] to [0, 2 J V], satisfying 0j(O) = and 6j(T) = 2- 7 7r. Since the Berry phase only depends 
on the path of the evolution of the Hamiltonian, the choice of the function 8(t) does not affect our result. Combining 
Eq. (fTB"]) and Eq. (|T5|) , the Hamiltonian of the j-th evolution is 



Hifii) =1- i(l+cos%)|^ )(Vo| - isin^-IViXV'il 

- 1(1 - 008^)1^X^1 + 1 sin^sK^I ! " i! 

7-|^)>(V^)l, 



where 

|^))=V^|6) + e-^V^|l), (17) 

and 9j G [0,2%]. Clearly, the ground state of H{9j) is \tp(6j)). We set the Hamiltonian of the whole system in the 
j-th adiabatic evolution as 

|0)(0| ® H{0j) + |1)(1| ®H(-6j). (18) 
Then the Berry phases of the computing system after the evolution are and — 7j, where 



i/ (M—Mde = ad9 = 2 J Tia. (19) 



Since the dynamic phase is the same in both parts, the final state is 

1 

3 



l^-(T)) = — (e^|0) + e -<7i|l))® |^ ). (20) 
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Hence, at the end of the evolution, the relative phase of the first qubit is Tj — 2jj — 2-k{2?o). In this way, we 
successfully encode the solution to the relative phases of a set of quantum states. 

In the algorithm, we do not use the phase estimation procedure in [13] to estimate a from the Berry phases, in 
order to avoid unnecessary entanglement. We use Kitaev's equivalent procedure instead [1 81 ] . As described above, 
in the j-th adiabatic evolution of the algorithm, we prepare a quantum state whose relative phase is 2ir(2 3 a), for 
j = 1, • • • , m. A measurement for the first qubit in |+) = ^(|0) + |1)) and |— ) = -^(|0) — |1)) basis gives the result 
|+) with probability 

p = cos 2 {^f) = cos 2 (2 J ira) . (21) 

Hence, if we apply the process several times, we can approximate the probability. More precisely, let q — r/R be the 
ratio between the number of r of results |+) and the number R of measurements. Then Chernoff's bound 

Prob{\p -q\>6)< 2e- 2&2R (22) 

shows that for a fixed 6, the error is smaller than e for only 0(log(l/e)) number of measurements. In the application, 
we only need an error that is smaller than n/8. Then, we obtain an estimation of 2 J a, modulo 1, with error 1/16. Let 

oo 

a = J2 2 ~ j <*j> (23) 

for a, S {0, 1}, and at = since a < 1/2. We also use .a± ■ ■ ■ a p to denote the binary fraction Xft-i %~^ a j- For j = 
1, • • • , to, we replace the known approximate value of 2 J a by rjj, the closet number from the set {0, 1/8, 2/8, • • • , 7/8}. 
Hence, we have 

|2 J a- 7ft 1 1 < 1/16+ 1/16 = 1/8. (24) 

Since if \y — 2a\i < S < 1/2, then \y' Q — a\i < 8/2 or \y[ — a|i < (5/2, where y' 1 ,y' 2 are the solutions to the equation 
2y' = j/(modl), we can start from 2 m a and increase the precision in the following way: Set r\ m — .a m a m +iot m +2 = Vm 
and proceed by iteration: 



a = I if\.0oj+iOj+2 - »ft|i < 1/4, , 25 . 
3 \\ if|.la J - +1 a J+2 -r? J -| 1 <l/4, 

for j = m — 1, • • • ,1. By a simple induction, .aia 2 ■ ■ ■ a m can estimate a with error less than 2 _m = e. 



III. RUNNING TIME OF THE ADIABATIC ALGORITHM 



In this section, we consider the accuracy of the evolutions and the running time of the algorithm. It is easy to see 
that under the Hamiltonian given in Eq. (|16p the actual state in the computing system, \ip.{6)), always stays in the 
subspace spanned by {|0), |1)}. Then H(9) can be rewritten as 



H(9) = 




in the subspace. Assume ui = 4| is constant, and u<l. Let |<^(£)) be the i-time state in the system which is initially 
in |^o)> an d evolving under H(ujt). By solving the Schrodinger equation ([1]), we attain 

\<p(t)) = e' l ^{Ae wlt + Be luJ2t )\0) + e^iCe" 1 " 11 + De~ W2t )\l), (27) 

where 



u> ± - uj) 2 + 4auj 

Ul,2 = 7, 



(28) 
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and 

A = (jz ^) 2 - "(1 ~ 3^) + (13- u)E ^ 



(1 


- lo) 2 + Aaui + (P — a- 


-lo)E 




a(l + lo) — aE 




(1- 


- uj) 2 + 4auo + ((3 - a - 


lo)E V 


(1- 


t lo) 2 - fi(l + 3lo) ~ (a 


+ lo)E 


(1 


- uj) 2 + iaoj + ((3 - a- 


-lo)E 









(29) 



D 



(1 - lo) 2 + Aollo + ((3 -a- uj)E 



where E — (1 — lo) 2 + Aauj, see the Appendix. On the other hand, denote the quantum state evolving under the 
Hamiltonian H(—6) = H(—Lot) by \tp'(t)), which can be obtained from \ip) by exchanging all to by —lo. 
Then, the final state of the j-th evolution is 



1 

V2 V 



\^{T)) = — (|0)|^(T)) + |1)|^(T))), (30) 



where T — 2 : >tt/lo, which is the time of the j-th evolution, for j = 1, • • • , m. As indicated before, we will measure the 
relative phase of the first qubit, and use it as an approximation of 2i:(2 :, a). Let \tp±) be a state in the span space of 
|0) and |1), which is orthogonal to \ip), from Eg . (f30|) . 



|r,,T» =-L(;0) + ( V '(T)y(T))il))^)~-l=(( V ' / / ! / C!l) 



Hence, with probability 



l + \( V (T)W(T))\ 2 
Ps = 2 ' ^ > 



the relative phase will be the argument of (ip(T)\(p' (T)) . It can be checked that 

(ip(T)\<p' (T)) = (AA' + DD')e ifil + {BB' + CC')e~ 1 ^ + (AB' + C'D)e if * 2 + (A'B + CD') e -^ 2 , (33) 

where A', B', C , D' are derived from A, B, C, D by exchanging to by — lo, and 

1 



M i = -( - lo) 2 + Aauo - + lo) 2 - 4aw)T, 



H2 = \W(^ - ^) 2 + 4au; + y/{\ + oj) 2 - 4aoo)T 
From the assumption lo <C 1, we have 

AA 1 + DD' = 1 - 3af3uo 2 + 0{lo 3 ), 
BB' + CC = -a/3L0 2 + 0(lo 3 ), 
AB' + C'D = 2a(3to 2 + 0(lo 3 ), 
A'B + CD' = 2a/3Lo 2 + 0(lo 3 ). 



(34) 



(35) 



Hence, 



and 



Moreover, 



1 -±^M^l>l-Sa^, (36) 



\arg({<p(T)\<p'(T))) - Mil < Sa(3co 2 . (37) 



Hi = 2auT + 2af3{(3 - a)uo 3 T + 0(lo 5 T) 
= 27r(2^a)(l + /3(/3 - a)co 2 + 0(lo 3 )), 



(38) 
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in the j-th evolution. Hence, the difference between the expected relative phase and 2ix(2 3 a) is 

A < 2tt(1 - Ps ) + \arg((<p+(T)\<p-(T))) - ,n| + |^i - 27r(2*a)| = 0(Vuj 2 ). (39) 

For the estimation of 2- J a modulo 1 with error less than 1/16, it is enough to make A < 27r/32, and then estimate 
the relative phase to its expected value within error 2n/32. To satisfy the first condition, in the j-th evolution we set 

u = Uj = 0(^=) (40) 



in Eq. [35] When estimating the relative phase, we can boost the success probability using repetitions. In our case, 
0(m — j) repetitions yield overall success probability greater than 1/2. Since the time of the j-th evolution is , 
the total running time is 

m 03 1 

T total =yo(-(m-j)) = 0((2" 1 ) 3 / 2 ) = 0((-f' 2 ). (41) 



Therefore, the adiabatic algorithm has an 0((i) 3 ^ 2 ) running time. As we know the optimal quantum algorithm 
has a running time O(^), our adiabatic algorithm is not as good as the quantum algorithm given in the circuit model. 
However, it is better than the best classical algorithm, which needs 0((|) 2 ) running time. 

IV. CONCLUSION 

In conclusion, we have proposed an adiabatic algorithm for the quantum counting problem. The key idea of the 
algorithm is to construct Berry phases which equal 2n(2 :, a), for j = 1, • • • , m, where a is the proportion of marked 
items in the database, m = log(-l). The algorithm has a running time of 0((i) 3 / 2 ), which beats the classical random 
algorithm. 

There are some special characteristics of our algorithm. Firstly, the solution of the problem is encoded in the phase 
difference of the final state, rather than the ground state of the final Hamiltonian. The final information is achieved 
by measuring relative phase, rather than the usual quantum measurement. Next, different from usual adiabatic 
algorithms, such as in [3j, we use more than one Hamiltonian as oracles to construct the evolution path. The idea of 
designing adiabatic algorithms based on geometric Berry phases could also be applied to solve other problems, such 
as searching an unordered database (2lj . 

We are grateful to Joseph F. Traub, Henryk Wozniakowski, Columbia University, Luming Duan, Yongjian Han, 
University of Michigan, and Mingsheng Ying, Tsinghua University, for their very helpful discussions and comments. 

V. APPENDIX 

In the Appendix, we will show the details of how to derive the quantum state under the Hamiltonian H (8) = H(u>t). 
The quantum system is initially in a state \ifj) = (y/]3, y/ct), then evolves under the Hamiltonian 

H $) = ( /^L-tf ~^ el9 \ . (42) 



Let \ f(t)} be the i-time quantum state in the system. Then the Schrodinger equation |T]) turns out to be, 

*>=^=(-^-- "^V- (43) 

Let \tp) = ( ) , then 



i^x = ax- y/^f3e luJt y, (44) 



ij t y = -^j3e-^ t x + f3y. (45) 
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From Eq|44l 



P V&P 



, dx 
~dt' 



dt 



lap 



\-iaux+{*-u)--i^\ 



Substituting Eqf46] into Eql45l 



Similarly, we can derive 



d 2 x dx 

—r + ill — UJj— + aujx = 0. 
dt z dt 



d 2 y 



dy 



dt2+l (l + .)--p.y = 0. 



The roots of Eq ATI and Eq[3H]are — § ± iu)\ % separately, where 



1 ± — oj) 2 + 4aw 



Then 



-i\t ( a iuji t 



< " 2 (Ae 
y = e"' = 



Be^ 2t ) 



From the initial state — \/]3\0) + \/a|l), we have 

A + B = V /^;C + L»- V / ^- 
Substituting Eqj5D] into EqHH it can be derived that 



D 
A 



B 

C 



Then, we attain that 



A = 



B = X 

c = i 



(/3 - a - w) + - uj) 2 + Aaui 

X^+y/P _ (1 - uj) 2 - a(l - 3a;) + (g - a^g 
1 + A 2 " 
Aa//3 — v 7 ^ 



D = X 



(1 - w) 2 + 4aw + (P - a - <jj)E 
a(l + to) — aE 
1 + A 2 ~~ (1 - uj) 2 + Aato + {fl - a - uj)E 
a - X^/P _ (1 + uj) 2 - P{\ + 3cj) - (a + 

(1 - w) 2 + 4aw + (/?-«- w)£ 

/?(l-aQ + /3£ 
' (1 - to) 2 + Aauo + (P-a- uj)E 



VP, 
y/P, 



1 + A 2 
AV^+ V/3 



1 + A 2 



(46) 



(47) 



(48) 



(49) 



(50) 



(51) 



(52) 



(53) 



where E = y(T~— a;) 2 + 4aw. On the other hand, denote the quantum state evolving under the Hamiltonian H(—9) 
H(—ujt) by \tp'(t)), and it can be derived from \tp) by exchanging all to by —uj. 
Then, the final state of the j-th evolution is 



|^(T)) = — doi^ + ii)!^™, 



(54) 



where T = 2- J 7r/o;, for j = 1, • • • , m. Let \<p±) be a state in the span space of |0) and |1), which is orthogonal to \cp), 
from Ea.([5i)). 



1 



1 



l^-(T)) = -=(|0) + (<p(T)\<S(T))\m<P) + ^=(<P±(TW(T)}\l)\v±(T)}. 



V2 



V2' 



(55) 
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Hence, with probability 



* - 1 + l( " (T f ,(r » |2 . (56) 



the relative phase will be the argument of (ip(T)\tp' (T)) . Since 

(<p(T)\(f'(T)) = (AA' + DD'y* 1 + (BB' + CC^e'^ 1 + (AB' + C D)e^ 2 + (A'B + CD^e'^ 2 , (57) 
where A', B' , C, D' are derived from A, B, C, D by exchanging uj by —uj, and 

Ml = 1( V(l - + Aauj - V(l + u) 2 - Aauj)T, 

2 (58) 
H2 = -(V(l - ^) 2 + 4aw + V(l + - 4a^)T 

From the assumption w<l, the relative phase can be estimated by {ip(T)\ip' (T)} , from AA' + DD' } BB 1 + CC, 
AB' + CD and A'B + CD'. 

Since the above 4 terms share the same denominator, we can first calculate the denominator F, then calculate their 
numerators separately. The denominator is 



F = [(1 - uj) 2 + Aauj + (1 - 2a - uj) yj (1 - uj) 2 + Aau] [(1 + uj) 2 - Aauj + (1 - 2a + uj) + uj) 2 - Aauj] 

= [(1 - uj) 2 + Aauj] ■ [(1 + uj) 2 - Aauj] + [(1 + uj) 2 - 4aw](l -2a- w)(l - (1 - 2a)uj + 2a0uj 2 ) 

+ [(1 - uj) 2 + Aauj](l - 2a + w)(l + (1 - 2a)uj + 2a0uj 2 ) r 

(59) 

+ [(1 - 2a) 2 - uj 2 ](1 - (1 - 2a)uj + 2a0uj 2 )(l + (1 - 2a)uo + 2a0u 2 ) + 0(uj 3 ) 

= 1 - 2(1 - 8a0)uj 2 + 2(1 - 2a) [1 - 2(1 - ha0)uj 2 ] + (1 - 2a) 2 - uj 2 - {I - 2a) 2 (l - 8a0)uj 2 + 0(uj 3 ) 
= A0 2 (1 - 2/3(1 - Aa)uj 2 ) + 0(uj 3 ). 

Then 

(AA' + DD')F 

= 0[(l - uj) 2 - a(l - 3w) + (1 - a - ^^(1 - u) 2 + 4aw] • [(1 + uj) 2 - a(l + 3cj) + (1 - a + uj) yj(l + uj) 2 - Aau)} 



+ a0 2 [(l -uj) + V '(1 - oj) 2 + Aauj] • [(1 + uj) + y/(l + uj) 2 - Aauj] 

= (3[(l - uj 2 ) 2 - 2a(l - huj 2 ) + a 2 (l - 9uj 2 ) + a0(l - uj 2 )] + 0[0 + /3(2(3 - 1)uj + (2a - 1)uj 2 ](1 - (1 - 2a)uj + 2a0uj 2 ) 
+ 0[0 - 0(20 - l)uj + (2a - 1)uj 2 ](1 + (1 - 2a)uj + 2a0uj 2 ) 

+ 0(0 - uj 2 )(2a - 1)uj 2 ](1 - (1 - 2a)uj + 2a0uj 2 )(l + (1 - 2o)uj + 2a0uj 2 ) + 0(uj 3 ) 
= A0 2 (1 - 0(2 - ba)uj 2 ) + O(uo 3 ). 

(60) 

Hence, 

+ dd> = !fJ!'5l'^l + Tl - 1 - 3^ 2 + o^ 3 ). (6i) 



4/3 2 (l - 20(1 - Aa)u 2 ) + 0(uj 3 ) 



(BB' + CC')F 



2 0(1 +uj- V (I - w) 2 + 4aw)(l - w - v / (l + ^) 2 -4aw) 



+ a[(l + w) 2 - + 3uj) - (a + uj) yj (1 - uj) 2 + Aau] ■ [(1 - uj) 2 - 0(1 - 3uj) -(a- uj)y/(l + uj) 2 - Aauj] 

= a[a + (9a0 + 2 - 2)uj 2 ] + [a + (1 - 2a)auj + (1 - 2a)w 2 ](l - (1 - 2a)uj + 2a0uj 2 ) (62) 

+ [a - (1 - 2a)auj + (1 - 2a)a> 2 ](l + (1 - 2a)uj + 2a0uj 2 ) 

+ a(a - uj 2 )(1 - (1 - 2a)uj + 2a0uj 2 ) • (1 + (1 - 2a)uj + 2a0uj 2 ) + 0(uj 3 ) 

= -Aa0 3 + 0(uj 3 ). 

Hence, 

BB ' + cc ' - m. -^f-l^V om - -°^ 2 + 0( " 3) - (63) 
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(AB' + C'D)F 

= a/3[(l - uj) 2 - a(l - 3oj) + {(3- uj)^ (I ~ uj) 2 + Aauo] • (1 - u - (I + uj) 2 - Aauj) 
+ a@[(l - lu) 2 - /3(1 - 3lu) + (a- + lu) 2 - Aauj] ■ (1 - w + ^ (I - uj) 2 + Aauj) 

= a/?(l - w)(l - w + 2lu 2 ) + a/3(l - (3 - 2/3)w + 2w 2 )(l - (1 - 2a)w + 2a/3w 2 ) (64) 

- a(3(l + (3 - 2/3)w + 2w 2 )(l + (1 - 2a)w + 2a(3uj 2 ) 

- (1 - 2w)(l - (1 - 2a)w + 2af3u?) • (1 + (1 - 2o)uj + 2af3uj 2 ) + 0(u 3 ) 
= 8af3 3 L0 3 + Q(uj 3 ). 



Hence, 



- , , 8aB 3 uj 3 + 0(uj 3 ) o , 

AB' + C'D = —— A . V 9 / — = 2aBuj 2 + (uj 3 ). 65 

AB 2 (1 - 2(3(1 -Aa)uj 2 )+0(uj z ) v ' v ' 

Since A'B + CD' = (AB' + CD)', 

A'B + CD' = 2aBuj 2 + 0(uj 3 ). (66) 

Hence, 

(tf(T)\tp'(T)) = (1 - 3a/3w 2 )e Vl - a^e^ 1 
+ AaBuj 2 cos(^ 2 ) + 0{uj 3 ). 

So 



(67) 



and 

|or fl ((p(T)b'(T))) - Mil < 8a/3^ 2 . (69) 
T = 0((i) 3/2 ) (70) 
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